Spin polarization control via magnetic barriers and spin-orbit effects 
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We investigate the spin-dependent transport properties of two-dimensional electron gas (2DEG) 
systems formed in diluted magnetic semiconductors and in the presence of Rashba spin-orbit inter- 
action in the framework of the scattering matrix approach. We focus on nanostructures consisting 
of realistic magnetic barriers produced by the deposition of ferromagnetic strips on heterostructures. 
We calculate spin-dependente conductance of such barrier systems and show that the magnetization 
pattern of the strips, the tunable spin-orbit coupling, and the enhanced Zeeman splitting have a 
strong effect on the conductance of the structure. We describe how these effects can be employed 
in the efficient control of spin polarization via the application of moderate fields. 
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Spin-orbit coupling in semiconductors intrinsically 
connects the spin of an electron to its momentum, ^ pro- 
viding a pathway for electrically initializing and manip- 
ulating electron spins for applications in spintronics^i^ 
and spin-based quantum information processing.^ This 
coupling can be regulated with quantum confinement in 
semiconductor heterostructures through band structure 
engineering, as well as by the application of external elec- 
tric fields, as in the celebrated spin field-eff'ect transis- 
tor proposed by Datta and Das.^ Using diluted magnetic 
semiconductors (DMS) in such systems provides an ad- 
ditional degree of control of the transport properties. In 
particular, when an external magnetic field is applied, 
the magnetic dopant spins align, giving rise to a strong 
exchange field that acts on the electron spin. This s- 
d exchange interaction between the electron spin in the 
conduction band and the localized magnetic ions induces 
a giant Zeeman splitting. 

In this communication we investigate the spin- 
dependent transport properties of two-dimensional elec- 
tron gas (2DEG) systems formed in diluted magnetic 
semiconductors and take into account the electric-field- 
dependent Rashba spin-orbit (SOI) interaction. We fo- 
cus our attention on nanostructures consisting of realis- 
tic magnetic barriers produced by the deposition of fer- 
romagnetic strips near the heterostructures;^ providing 
a relatively strong inhomogeneous magnetic field on the 
2DEG. We show how the conductance of the 2DEG de- 
pends strongly on the magnetization pattern of the strips, 
as well as on geometry and externally applied electric 
fields. We demonstrate that significant spin polarization 
(exceeding 50%) can be obtained at low temperatures for 
ferromagnetic strips of typical dimensions and magneti- 
zation. 

Figure [1] illustrates the type of magnetic barrier which 
can be created by the deposition of a ferromagnetic strip 
on the surface of a heterostructures. In this case the 
magnetization is assumed perpendicular to the 2DEG lo- 
cated at a distance zq below the surface. For this strip 



of width d and thickness h, the magnetic field along the 
z-axis is given by^ 



B = Bo[f{x + d/2)-f{x-d/2)], 



(1) 



where Bq = Moh/d, f{x) = 2xd/{x'^ + Zq), and Mq 
is the magnetization of the ferromagnetic strip. The 
in-plane field component can be assumed negligible for 
this configuration)^ We should mention that deposition 
of dysprosium strips result in magnetic fields estimated 
to be as large as IT in the vicinity of the 2DEG;^ from 
this we take B^ = O.IT to be a realistic value. 

The Hamiltonian of the 2DEG (in the xy plane) is 
realized in the lowest subband of the semiconductor 
heterostructure with effective mass m*, and includes a 
Rashba SOI due to z confinement characterized by the 
coupling constant a^, Zeeman term with electronic g fac- 
tor, and the s-d exchange interaction: 



B{x) + J2 «(^) • SiR^)Sir ~ R,), (2) 



where 11^ denotes the kinetic momentum, and cr^ the 
Pauli spin matrices. S is the spin of the localized 3d^ 
electrons of Mn^+ ions with S = 5/2, and s is the elec- 
tron spin in the 2DEG. We assume that the magnetic ions 
are distributed homogeneously in the DMS, and that the 
extended nature of the electronic wave function spans 
a large number of magnetic ions, allowing the use of a 
molecular-field approximation to replace the magnetic- 
ion spin operator Si with the thermal and spatial av- 
erage (Sz) along the external magnetic field direction. 
This approach has been proven to be suitable in previous 
studieSi^ii^ Correspondingly, the exchange interaction in 
Eq. ([2]) can be written as a Zeeman-like term, Jsd{Sz)o'z, 
where (S,) = {5/2)Bj{SgfiBB/kB{T + To)), with Bj{x) 
as the Brillouin function, Tq accounts for the reduced 
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FIG. 1: (color online) Top diagram illustrates magnetic strip 
with perpendicular magnetization over a 2DEG layer. The 
magnetization of the strip results in the magnetic barrier and 
vector potential shown in middle panel. Bottom panel shows 
effective potential Ve for different q values (Eq. Structural 
parameters chosen as d = 1 and zq = 0.1, in units of the 
magnetic length Lq ~ yfi/sSo =8.1 nm for Bq = O.IT. 



single-ion contribution due to the antiferromagnetic Mn- 
Mn coupling, and ks is the Boltzmann constant. The 
parameters used in the calculation are Jsd = —0.22a; eV, 
X = 0.014, To = 40K,W and T = IK. For convenience 
we express all quantities in dimensionless units, using 
Lo = eBo/m*, and magnetic length Lq = ^JK/cBq. For 
Gai_a;Mn2;As, m* = 0.067me, one gets in Lq — 8.1 nm 
and huj = 0.17 meV, for Bq — 0.1 T as above. 

The two-dimensional Schrodinger equation H'^ = E'i 



Fig. [H it is difficult to solve Eq. (|4]) analytically. How- 
ever, a numerical approach which divides the region 
of the barrier [a;i,a;/] into 7V(^ 1) segments of width 
L = {xf — Xi)/N is possible.^ For small L (large N) the 
vector potential in each of the segments can be treated as 
constant, so that in each jth segment one has solutions 



liven by i}j{x)i 



where 



X±{kj), with spinors 



xl 



(5) 



and 



A. 



-at[q + A{xj) + ik,] 
v/2A±(A± - Z,) 

v/2A±(A±-^,) ' 
J:AS.) + 19*B{x,), 



±{Z] + \al[q + A{x,) + ik,]\^} 



1/2 



(6) 



(7) 



where kj is the solution of the equation [k'j+{q+A{xj))'^~ 
2EY = Aj_. We consider the transmission of electrons 
with initial wave vectors A;"*" and k~ corresponding to the 
region x ^ Xi, which are the solutions of the problem 
with Jsd{Sz) — 0, and = away from the magnetic 
barrier region. It is clear that the transmission probabil- 
ities for electrons of spin eigenstates and "— " in ([5]) 
are different due to the SOI, as well as by the enhanced 
and inhomogeneous Zeeman splitting in Since the 
wave vector in the y direction is fixed, on both inci- 
dent and outgoing regions of the magnetic barrier, the 
wave function can be written as i^ix) = e*'^^^x+(fc"'") -I- 

X~{-k-) 
X-{k~) for 



for X < Xi\ 



\-k+) + e**="^x-(fc") + 
and as t'^e^'^'^^x+ik'^) + 



-ik 



X > Xf. Here and are the transmission and re- 
flection coefficients for each spin state x±- Taking into 
account the boundary conditions, which require that the 
wave function be continuous at the interfaces, a system 



has solutions of the form *(x, y) = e^^y^^ix), where E oflinear equations for t± and r± can be derived. Once t 



is the total energy of the electron, and q is the electron 
wave vector in y direction. 4^{x) satisfies 



dx^ 



Ve{x,az) + al[axiq + A{x)) (3) 



'<Jy{~t-^)]\^/j{x)^2E^/j{x) 



Ve{x, a,) = [q + A{x)]^ - [J:,{S,) + -9*B{x)]a, (4) 

where a* = 2az/LohLL>, J*^ = 2Jsd/f^, 9* = 9^iBBa/hhJ, 
and Ve is the effective dynamical potential. We use 
the Landau-gauge, A{x) = (0,A(a;),0), and B{x) = 
dA{x)/dx. 

For a magnetic barrier structure as that shown in 



and r are known, it is straightforward to obtain the spin 
dependent transmission coefficients Tss> {E, q) , where s' 
(s) is the incident (outgoing) spin, with s, s' =t or |. 
Note that for incident electrons with spin s' =t there 
are two kinds of transmission coefficients, Tji, and T|f , 
as the SOI produces precession of the electron spin as it 
propagates. A similar effect occurs for incident electrons 
with spin s' —[, although the presence of the magnetic 
field breaks time-reversal invariance which results in a 
spin- filtering effect, as we will see below. 

In the ballistic regime the spin-dependent conductance 



GsAT) 



9sAE) - 



d fiE,T) 
dT 



dE, 



(8) 



3 



^ 0.4 
0.2 

_aj 

Ml 

g 

1-0.2 
u 

^-0.4 

^ 0.4 

0.2 

_aj 

Ml 

^ 

1-0.2 
u 

-S-0.4 



I 



(b) Til 



0.8 



0.6 



0.4 



0.2 



20 



40 60 
Fermi Energy {huj) 



80 



100 



FIG. 2: (color online) The transmission coefficients (a) T||, 
(b) as function of energy and incident angle 9 for electrons 
tunneling through the magnetic barrier with perpendicular 
magnetization, d = 2Lo and SOI coupling Oz = 0.35 x 10~^^ 
eV-m. 



where 



gss' (T) = Go 



/ T,s'{E, V2E sin{e)) cosiO) d0 , (9) 

J-Tl/2 



9 is the angle of incidence with respect to the x-direction 
(see top panel Fig. [1]), f{E,T) is the Fermi-Dirac distri- 
bution function, and Go = e'^m*VFLy/h^ — [e^ /h)kFLy, 
where Ly is the length of the structure in y direction and 
vp = kkp is the Fermi velocity. 

To evaluate the electron spin polarization effect in the 
tunneling process, we define the net conductance polar- 
ization as 
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In Fig. [5] we show typical transmission coefficients, as 
function of incident energy and angle [q — \J2E sin 6) , 
for an incoming electron with spin |. The upper panel 
shows the transmission coefhcient for collected electrons 
with the same polarization, T^^, while the lower is for 
outgoing electrons with opposite spins, T||. Here we use 
a typical SOI strength, = 0.35 x 10^^^ eV-m,'-^ which 
corresponds to a* = 0.5 . Notice that for low incident 
energy and/or with a large incident angle 9 the transmis- 
sion coefficients are very small, as the effective barrier 
presented by the magnetic strip is quite large (see bot- 
tom panel of Fig. [T|). At lower incident angle and higher 
energy, the transmission exhibits a series of resonances 
near normal incidence (9 « 0), as the magnetic vector 
potential results in an effective scattering double barrier 
potential (14 ~ A^[x) in Eq. 0]) with concomitant res- 
onances. This structure is also seen in the spin mixing 
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FIG. 3: (color online) (a) Spin-dependent conductance (in 
units of Go) and (b) polarization as function of Fermi energy 
(in units of hui). Structure parameters as in Fig. [2] 



results, T||, as well as in the corresponding spin | in- 
cidence. Til and T^i (not shown), although with some 
subtle differences due to the asymmetry introduced by 
the field. 

The calculated conductance and polarization compo- 
nents for the same parameters as in Fig. [2] are shown 
in Fig. [21 vs. Fermi energy. The SOI naturally causes 
the spin of the electron to precess when the electron 
propagates through the gated region. This is evident in 
Fig. [SKa) as the spin up conductance G^i is quite differ- 
ent from he spin down conductance, G||. The inhomo- 
geneous magnetic field also contributes significantly to 
change the tunneling probability of the electron through 
the magnetic gate barrier. G|| exhibits peaks shifted 
in general towards lower energy values, in contrast with 
those in Gn, a behavior persisting at low energies, where 
the two curves show clearly split resonances. In addi- 
tion, the spin mixing probability of tunneling electrons 
depends greatly on the incident spin for a given mag- 
netic barrier, as illustrated by the curve for G|| being 
slightly different to G|| (see especially the insert in Fig. 
Ha)). 

The conductance polarization plotted in Fig. ^h) 
shows that this structure results in polarization over 50% 
in the first resonant energy and close to 70% at the sec- 
ond and higher resonances. Therefore, for realistic values 
of strip magnetization and spin-obit coupling the system 
can generate substantial spin polarized currents, even as 
the injection is unpolarized. Notice that the resonant 
peaks created by the effective dynamic potential are re- 
sponsible for the large polarization values, as the latter 
decreases rapidly for higher Fermi energies. 

To provide a more complete account of the role of SOI 
in this polarization effect we show in Fig. 2] the calcu- 
lated conductance and polarization as functions of the 
SOI strength a* at two different Fermi energies. [Notice 
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FIG. 4: (color online) Spin conductance and polarization as 
function of for different Fermi energies. Notice polarization 
exceeds 50% in panel (d) and reverses in a small a* interval. 



that a can be varied in experiments by the application 
of a gate voltage applied to the metallic strip itself.] Fig- 
ure |4] (a) and (b) show results for Ep = 63?iw. At this 
relatively large energy the spin dependent conductances 
show regular oscillations with a*, similar to those ex- 
pected from the Datta-Das devicci^ In this region of SOI 
the polarization starts near zero (due to Zeeman effect) 
and has a local maximum at a* = 0.45 of no larger than 
1%, becoming slightly negative afterwards. Figure 2] (c) 
and (d) show results for Ep = 27.8hLLj, near one of the 
transmission resonances in in Fig. [3] It can be seen 
that increasing a results also in the spin conductances 
having oscillations. These features reflect the spin pre- 
cession effects induced by the SOI, further enhanced by 
the resonances of the dynamical potential, as discussed 
above. Notice that in this case, the conductance polar- 
ization can however exceed 50% and be fully reversed 
when a* changes by ~ 0.23. 

As the strong spin polarization of the conductance is 
associated with the presence of transmission resonances 
through this combination of magnetic and electric field 
barrier, we also explore its dependence on strip width- 
which determines the width of the effective double barrier 
potential 14 . In Fig. [5] we show the spin polarization vs. 



Fermi energy for different magnitudes of the strip width 
d {d = Lq to 4Lo)- One can see that increasing d results 
in a larger width of the double barrier potential well pro- 
duced by the effective potential, which in turn results in 
more resonant peaks in the conductance as the Fermi en- 
ergy increases. In all cases, notice that the polarization 
exceeds 50% to 70%. 

We have shown that it is possible to achieve a sub- 
stantial spin polarization in a realistic situation where a 
magnetized strip is placed in close proximity to a 2DEG. 
Our calculation incorporates the role of the giant Zee- 
man effect in diluted magnetic semiconductors, and this, 
together with the tunable SOI coupling, can be appro- 
priately used to enhance or suppress the polarization. 
We should also notice that as the temperature changes. 
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FIG. 5: (color online) Spin polarization as function of Fermi 
energy at a given value of a* — 0.5 for several width values, 
d, of the magnetic strip. Curves offset for clarity. 



the effective field can be modulated, with the consequent 
changes on the spin polarization characteristics of a given 
structure. 

The authors thank the support of CMSS and BNNT 
programs at Ohio University, as well as NSF grants 
0336431 and 0710581. JVB acknowledges support of the 
A. von Humboldt Foundation. We thank helpful conver- 
sations with D. Csontos and M. Zarea. 



^ R. Winkler, Spin-Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems, Springer Tracts 
in Modern Physics (Springer, Berlin, 2003) Vol. 191. 

^ S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. 
Daughton, S. von Molnar, M. L. Roukes, A. Y. Chtchelka- 
nova and D. M. Treger, Science 294, 1488 (2001). 

^ I. Zutic, J. Fabian and S. Das Sarma, Rev. Mod. Phys. 76, 
323 (2004). 

^ Semiconductor Spintronics and Quantum Computation, D. 
D. Awschalom, D. Loss, and N. Samarth, eds. (Springer, 
Berlin, 2002). 

^ S. Datta and B. Das, Appl. Phys. Lett. 56, 665 (1990). 



^ A. Matulis, F. M. Peeters and P. Vasilopoulos, Phy. Rev. 
Lett. 72, 1518 (1994). 

P. D. Ye, D. Weiss, R. R. Gerhardts, M. Seeger, K. von 
Klitzing, K. Eberl and H. Nickel, Phy. Rev. Lett. 74, 3013 
(1995). 

J. Q.You, L. Zhang, P. K. Ghosh, Phys. Rev. B 52,17243 
(1995). 

^ A. Lemaitre, C. Testelin, C. Rigaux, T. Wojtowicz and G. 
Karczewski, Phys. Rev. B 62, 5059 (2000). 
J. C. Egues, Phys. Rev. Lett. 80, 4578 (1998). 
R. Lang, A. Winter, H. Pascher, H. Krenn, X. Liu, and J. 
K. Furdyna, Phys. Rev. B 72, 024430 (2005). 



Destefani, S. E. Ulloa, and G. E. Marques, Phys. Rev. B 70, 205315 (2004). 



